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I. INTRODUCTION 

The preferred model of cosmology today, the ACDM 
model, provides an excellent fit to cosmological observa- 
tions, but at a substantial cost: according to this model, 
about 95% of the universe is either invisible or unde- 
tectable, or possibly both [1]. This fact provides a strong 
incentive to seek alternative explanations that can ac- 
count for cosmological observations without resorting to 
dark matter or Einstein's cosmological constant. 

For gravitational theories designed to challenge the 
ACDM model, the bar is set increasingly higher by recent 
discoveries. Not only do such theories have to explain 
successfully the velocity dispersions, rotational curves, 
and gravitational lensing of galaxies and galaxy clus- 
ters, the theories must also be in accord with cosmolog- 
ical observations, notably the acoustic power spectrum 
of the cosmic microwave background (CMB), the matter 
power spectrum of galaxies, and the recent observation 
of the luminosity-distance relationship of high-z super- 
novae, which is seen as evidence for "dark energy" . 

Modified Gravity (MOG) ^2] has been used success- 
fully to account for galaxy cluster masses the rota- 
tion curves of galaxies [1] , velocity dispersions of satellite 
galaxies 0], and globular clusters [1]. It was also used 
to offer an explanation for the Bullet Cluster Q without 
resorting to cold dark matter. 

Remarkably, MOG also meets the challenge posed by 
cosmological observations. In this paper, it is demon- 
strated that MOG produces an acoustic power spectrum, 
a matter power spectrum, and a luminosity-distance re- 
lationship that are in good agreement with observations, 
and require no dark matter nor Einstein's cosmological 
constant. 

In the arguments presented here, we rely on simplified 
analytical calculations. We are not advocating these as 
substitutes for an accurate numerical analysis. However, 
a thorough numerical analysis requires significant time 
and resources; before these are committed, it is useful to 
be able to demonstrate if a theory is viable, and if the 
additional effort is warranted. 



In the next section, we review the key features of MOG. 
This is followed by sections presenting detailed calcula- 
tions for the luminosity-distance relationship of high-z 
supernovae, the acoustic power spectrum of the CMB, 
and the galaxy power spectrum. A concluding section 
summarizes our results and maps out future steps. 

II. MODIFIED GRAVITY THEORY 

Modified Gravity (MOG) is a fully relativistic theory of 
gravitation that is derived from a relativistic action prin- 
ciple [2] involving scalar, tensor, and vector fields. MOG 
has evolved as a result of investigations of Nonsymmet- 
ric Gravity Theory (NGT) [H, and most recently, it has 
taken the form of Scalar- Tensor- Vector Gravity (STVG) 
In the weak field approximation, STVG, NGT, and 
Metric-Skew- Tensor Gravity (MSTG) 9] produce similar 
results. 



A. Scalar-Tensor- Vector Gravity 

Our modified gravity theory is based on postulating 
the existence of a massive vector field, 0^ . The choice of 
a massive vector field is motivated by our desire to in- 
troduce a repulsive modification of the law of gravitation 
at short range. The vector field is coupled universally 
to matter. The theory, therefore, has three constants: 
in addition to the gravitational constant G, we must also 
consider the coupling constant lo that determines the cou- 
pling strength between the cj)^ field and matter, and a 
further constant fi that arises as a result of considering a 
vector field of non-zero mass, and controls the coupling 
range. The theory promotes G, /i, and uj to scalar fields, 
hence they are allowed to run, resulting in the following 
action @,[3: 

S = Sg + Sc/, + Ss + Sm, (1) 

where 
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where S'm is the "matter" action, Bni, = dn, 



d^(j)^, while Vc/,{(j)), Vg(G), V^{uj), and denote the self- 



interaction potentials associated with the vector field and the three scalar fields. The symbol is used to denote 
covariant differentiation with respect to the metric g^"^ , while the symbols R, A, and g represent the Ricci-scalar, the 
cosmological constant, and the determinant of the metric tensor, respectively. We define the Ricci tensor as 
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Our units arc such that the speed of light, c = f ; we use the metric signature (+, — , — , — ). 
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B. The field equations 



Using the FLRW line element ds^ = df - a(t)2[(l - kr^y^dr^ + r^d^^\ with dfi 
Friedmann equations in the form [lol |: 
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where H = a/a, and the dot denotes differentiation with respect to t, i.e., y — dy/dt. These equations are supplemented 
by an additional set of three equations for the scalar fields G, /i, and lj: 
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On the right-hand side of the Friedmann equations ^ 
and ([7]), in addition to terms describing ordinary matter 
and energy (characterized by p and p; herein, we are only 
considering the case of pressureless dust, such that p = 
and w — p/ p — 0), we see several additional terms, each 
with its distinct equation of state. We can rewrite the 
Friedmann equations as 
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The equations of state associated with the kinetic term 
Pk, potential term pv , cosmological term pA, and Brans- 
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Dicke term pQ are, respectively, 
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In particular, we note that the term associated with the 
potentials V^, Vg, V^, and has the equation of state 
wy = — 1, indicating that these potentials can play the 
role of dark energy, even if we set Einstein's cosmological 
constant, A = 0. 

Equations ([S HIP)) can be solved numerically, if suitable 
initial conditions are first established. We choose to use 
the following initial values: 

to = 13.7 X 10^ years, (21) 

ao = do, (22) 

Go = 6Gjv, (23) 

^0 - (24) 

LOO = l/%/T2, (25) 

00 = 0, (26) 

do = Hoao = 72ao km/Mpc s, (27) 

Go^fio=djo = 0, (28) 

Vg = 0.07659537Gg/tg, (29) 

V,, - = = 0, (30) 

A = 0, (31) 

fc = 0, (32) 

where Gn is Newton's gravitational constant. Most of 
these choices are self-explanatory, and consistent with 
standard cosmology. Exceptions are the scalar field val- 
ues. Our choices are motivated by the solution of the 
MOG field equations in the spherically symmetric case 
(for details, see fl^), yielding a Yukawa- like modifica- 
tion of gravity with /i as the range parameter. We set fi 
to the inverse of the scale of the universe. 

The value of Go is chosen to yield a model that requires 
no substantial amounts of nonbaryonic matter. With G 
at 6 times the Newtonian constant, baryonic matter alone 
amounts to ^ 30% of the critical density of the universe, 
hence no nonbaryonic dark matter is needed to explain 
a deficit in matter density. The value of 6Gn is also 
remarkable for another reason: taking the point particle 
solution from the next section, an effective gravitational 
constant of GeS — QGn at the Yukawa distance r = 
fi~^ yields a ~ 19 and an effective gravitational constant 
of Gcff — 20GAr at infinity. That is, on superhorizon 
scales our solution is consistent with an Einstein-de Sitter 
cosmology with no dark matter or dark energy. (This 
would imply a vanishing Vg on superhorizon scales.) 

With these initial values, the exact numerical solution 
of ([5 HTU)) yields a bouncing cosmology (Figure [T|) : the 
universe contracts up until approximately 0.9 times the 
Hubble age before the present epoch, at which time the 




FIG. 1: The MOG "bouncing" cosmology. The horizontal 
axis represents time, measured in Hubble units of Hq^. The 
solid (black) line is a/ao, the scale factor normalized to the 
present epoch. The dashed (green) line is G/Gq. The inset 
shows details of the bounce, demonstrating that a smooth 
bounce occurs even as the matter density of the universe is 
more than 10^* times its present value. 




FIG. 2: The effective equation of state w^s (solid green line) 
and the deceleration parameter q (dashed blue line) in MOG 
cosmology. Inset shows the behavior of the deceleration pa- 
rameter in the vicinity of the bounce. 



contraction halts and expansion begins. At this time, the 
value of G changes sharply, albeit smoothly; meanwhile, 
the scalar fields /i and u remain constant. 

This type of bouncing cosmology is known from the 
literature: Brans-Dickey theory |ll| pr oduces a smooth 
bounce if its parameter ojbd < ~6 |12l |. In our case, as 
the /i and lu fields are constant, our cosmology is similar 
to that of Brans-Dicke theory with wbd — — Stt. 

One feature not normally present in Brans-Dicke the- 
ory is the potential field Vg- Choosing an appropriate 
constant value for this potential allows us to fine-tune 
our cosmology, making the bounce sharper, yet keeping 
it smooth even as the universe achieves densities in excess 
of 10^^ times the present density. 

Furthermore, introducing a mixed equation of state for 
ordinary matter and radiation does not qualitatively alter 
this result: the essential features of the bounce as well as 
the time scales remain the same. 

To investigate further the features of the MOG cosmol- 
ogy, we define two additional parameters: the effective 
equation of state WcS and the deceleration parameter q. 

The effective equation of state is obtained by combin- 
ing all terms on the right-hand side of the Friedmann 



4 



equations ([6]) and ([7]) into an effective density pcS and ef- 
fective pressure p^g — WcsPcS, such that the Friedmann 
equations now read: 
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The deceleration parameter is defined in the usual way: 
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These two dimensionless parameters are shown in Fig- 
ure [H We note that the effective equation of state always 
remains confined between — 1 ^ Wcff ^ 0, except at the 
bounce. In the absence of curvature (fc — 0) and a cos- 
mological constant (A = 0), the bounce can only occur if 
d = 0, but at the same time, a must not be zero. This 
means that Poff is non-zero even as poff is zero, as can 
been from ([55)) and ([M)) . In other words, at the time of 
the bounce, as q goes to negative infinity, so does w^s- 

In the past, we have considered a bouncing cosmology 
in the MOG context [l^. We considered the possibility 
that at the time of the bounce, the universe was in a state 
of minimal entropy, which would define a thermodynam- 
ical arrow of time that points in the +t direction after 
the bounce, but in the —t direction prior to the bounce. 



C. Point particles in a spherically symmetric field 

For a point particle moving in the spherically sym- 
metric field of a gravitating source, a particularly simple 
solution for the acceleration is obtained Il4i: 
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where M is the source mass, while a determines the 
strength of the "fifth force" interaction, and /i controls 
its range. In prior work, a and /i were considered free 
parameters that were fitted to data. Our recent work [lo| 
allows us to determine a and [i as functions of the source 
mass M: 
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This solution can be seen to satisfy the field equations 
in the spherically symmetric case either numeric ally or by 
deriving an approximate solution analytically The 
numerical values for D and E are determined by match- 
ing the result against galaxy rotation curves [lOj : 
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(39) 
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The value of Goo — ^OGn is set to ensure that at the 
horizon distance, the effective strength of gravity is about 
6 times Gat, eliminating the need for cold dark matter 
in cosmological calculations, as described in the previous 
section. 



D. The MOG Poisson Equation 

The acceleration law ([55)1 is associated with the poten- 
tial, 
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is the Newtonian gravitational potential with Goo = (1" 
a) Gat as the gravitational constant, and 
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is the Yukawa-potential. These potentials are associ- 
ated with the corresponding Poisson and inhomogeneous 
Helmholtz equations, which are given by Q: 



V2$jv(r) = 4^GooP(r), 



(V2-A.2)$y(r) = 
Full solutions to these potentials are given by 
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These solutions can be verified against Eqs. and 
(|43p by applying the delta function point source density 
p{y) = M5^{r). 

Combining Eq. (gT]) with Eqs. (gl]) and (gS]) yields 



= 47rGArp(r) -h/i2$y(r) 



(48) 
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containing, in addition to the usual Newtonian term, a 
nonlocal source term on the right-hand side. 



III. MOG AND HIGH-Z SUPERNOVAE 

Type la supernovae are excellent standard candles for 
astronomy. The physics of these massive explosions is 
believed to be well understood, and it determines their 
peak luminosity. The distance to these events can be 
calculated from their redshift; knowing the distance and 
their absolute luminosity, their apparent luminosity can 
be computed. 
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FIG. 3: Type la supernova luminosity-redshift data [3| and 
the MOG/ACDM predictions. No astrophysical dimming was 
applied. The horizontal axis corresponds to the q = empty 
universe. The MOG result is represented by a thick (blue) 
line. Dashed (red) line is a matter-dominated Einstein de- 
Sitter universe with = 1, g = 0.5. Thin (black) line is the 
ACDM prediction. 



A. The luminosity-distance relationship 



and Vuj. While we have no a priori reason to single 
out Vc, the fact that it offers a means to "fine tune" 
the MOG bouncing cosmology is a major point favoring 
this choice. However, it must be said that interesting 
cosmologies can also be generated using the other self- 
interaction potentials. Conversely, the effects of Vg in 
the Friedmann equations are not unique to MOG, but 
also present in any other theory that contains a Brans- 
Dicke scalar field with negative wbd and a self-interaction 
potential. 

At the very least, our analysis shows that MOG is con- 
sistent with supernova observations. Uncertainties in the 
data, and questions concerning corrections such as as- 
trophysical dimming make it difficult to present a more 
definitive statement on the basis of these observations 
alone. 



The difference fi{z) between the absolute and appar- 
ent luminosity of a distant object can be calculated as a 
function of the redshift z, Hubble constant H, and cosmic 
deceleration q = — aa/a^ as [T5| : 

ti{z) = 25-5 logio H + 51ogio(cz) + 1.086(1 - q)z + ... 

(49) 

Even after corrections for astrophysical dimming are 
applied, supernova observations are not consistent with 
the deceleration parameter q — 0.5 value of a matter- 
dominated Einstein-de Sitter universe. 

The deceleration parameter can be expressed from the 
Friedmann equations ([55]) and as 



(? = -(1 + 3Weff) ( 1 + ^) - + ^cff);|2 ' (50) 



In the absence of exotic forms of matter, i.e., assuming 
that Weft > 0, and in the absence of curvature {k = 0), 
only a positive cosmological constant, A > 0, can reduce 
the value of 9 to a value that is consistent with supernova 
observations. 

In the case of MOG, however, the effective equation 
of state is WcS < 0, and descends below —1/3 when the 
universe is about two thirds its present age. This allows q 
to be negative, resulting in an accelerating universe even 
though A = 0. 



B. Discussion 

Consistency between its predictions and supernova ob- 
servations has been a major point in favor of acceptance 
of the ACDM model. However, the model leaves the 
question of the origin of the A term open: it may be 
Einstein's cosmological constant, it may be a potential 
associated with a scalar field (quintessence) , or it may be 
something else altogether, so long as it yields an equation 
of state w < —1/3. 

MOG provides just such a term in the Friedmann equa- 
tions in the form of the self-interaction potentials V^,, Vg, 



IV. MOG AND THE CMB 

The cosmic microwave background (CMB) is highly 
isotropic, showing only small temperature fluctuations 
as a function of sky direction. These fluctuations are not 
uniformly random; they show a distinct dependence on 
angular size, as has been demonstrated by the measure- 
ments of the Boomerang experiment (iTj and the Wilkin- 
son Microwave Anisotropy Probe (WMAP) fll|. 

The angular power spectrum of the CMB can be cal- 
culated in a variety of ways. The preferred method is to 
use numerical software, such as CMBFAST [l^. Unfortu- 
nately, such software packages cannot easily be adapted 
for use with MOG. Instead, at the present time we opt 
to use the excellent semi-analytical approximation de- 
veloped by Mukhanov [l^. While not as accurate as 
numerical software, it lends itself more easily to nontriv- 
ial modifications, as the physics remains evident in the 
equations. 



A. Semi-analytical estimation of CMB anisotropics 

Mukhanov [l^ calculates the correlation function C{1), 
where I is the multipole number, of the acoustic power 
spectrum of the CMB using the solution 



[C(0]iow; 9 ^^ + ''^' 



(51) 



where I ^ 1, O denotes the oscillating part of the spec- 
trum, while the non-oscillating part is written as the sum 
of three parts: 



N = Ni + N2 + N3. 
These, in turn, are expressed as 



A^i = 0.063C' 



[P-0.22{l/lfr^-2.6r n/^^y. 



H-0.65(V^/)i-4 



(52) 



, (53) 
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1.4 



^ 0.033 [P-Q.b{l/hr^^ + 2.2f _(;/,^). 
The oscillating part of the spectrum is written as 



, (54) 
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where 
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^{^pl + \)] , (56) 



Ao = 0.14 



(1 + 0^/^ 



(57) 



(58) 



The parameters that occur in these expressions are as 
follows. First, the baryon density parameter: 



e-i7(i)6/i?5), 



(59) 



where fif, ~ 0.035 is the baryon content of the uni- 
verse at present relative to the critical density, and 
/175 = H/{75 km/s/Mpc). The growth term of the trans- 
fer function is represented by 



P = In 



-0.09 



200v/aji 



(60) 



where flm — 0.3 is the total matter content (baryonic 
matter, neutrinos, and cold dark matter). The free- 
streaming and Silk damping scales are determined, re- 
spectively, by 

Ij: = 1300 [1 + 7.8 X 10"2 (nmh'^r,y^] ^'"^ n"J^, (61) 
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(62) 
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Lastly, the location of the acoustic peaks is determined 
by the parameter^ 



p = 0.015(1 + 0.13O"^(f^m/i?5^) 



3.1\0.16 



(63) 



^ Note that we slightly adjusted the coefficients of 11611 and (163 l l. 
which improved the fit noticeably, while remaining fully consis- 
tent with Mukhanov's derivation. 
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FIG. 4; MOG and the acoustic power spectrum. Calculated 
using Qm ~ 0.3, fit — 0.035, Ho — 71 km/s/Mpc. Also shown 
are the raw WMAP 3-year data set (light blue), binned aver- 
ages with horizontal and vertical error bars provided by the 
WMAP project (red), and data from the Boomerang experi- 
ment (green). 



B. The MOG CMB spectrum 

The semi-analytical approximation presented in the 
previous section can be adapted to the MOG case by 
making two important observations. 

First, in all expressions involving the value of 
Mukhanov's flm (which includes contributions from 
baryonic matter and cold dark matter using Newton's 
gravitational constant), we need to use — 0.3 (which 
includes baryonic matter only, using the running value 
of the gravitational constant, G^s — 6Gjv). Second, we 
notice that the value of Qb in (|59p does not depend on the 
effective value of the gravitational constant, as this value 
is a function of the speed of sound, which depends on 
the (baryonic) matter density, regardless of gravitation. 
In other words, ~ 0.035 is calculated using Newton's 
gravitational constant. 

After we modify Mukhanov's semi-analytical formu- 
lation by taking these considerations into account, we 
obtain the fit to the acoustic power spectrum shown in 
Figure SI 



C. Discussion 

As Figure 0] demonstrates, to the extent that 
Mukhanov's formulation is applicable to MOG, the 
theory achieves excellent agreement with the observed 
acoustic power spectrum. We wish to emphasize that 
this result was obtained without fine-tuning. The MOG 
constant /i was assumed to be equal to the inverse of the 
radius of the visible universe. Thereafter, the value of a 
is fixed if we wish to ensure JIm — 0.3. This was sufficient 
to achieve good agreement with the data. 
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V. MOG AND THE MATTER POWER 
SPECTRUM 

The distribution of mass in the universe is not uni- 
form. Due to gravitational self-attraction, matter tends 
to "clump" into ever denser concentrations, leaving large 
voids in between. In the early universe, this process is 
counteracted by pressure. The process is further compli- 
cated by the fact that in the early universe, the energy 
density of radiation was comparable to that of matter. 



Subtracting the original set of equations from the new 
set, using l/{p+5p) = {p-5p)/[p^-{5pf] = l/p-Sp/p^, 
and eliminating second-order terms, we obtain 

^ + V • (5pv + p5w) = 0, (66a) 



— ^-H(v-V)(5v+((5v-V)v = -^Vp--V5p + (5g, (66b) 
at p^ p 



A. Density fluctuations in Newtonian gravity 



V • (5g -47rG'(5/rj. 



(66c) 



To first order, this process can be investigated using 
perturbation theory. Taking an arbitrary initial distribu- 
tion, one can proceed to introduce small perturbations in 
the density, velocity, and acceleration fields. These lead 
to a second-order differential equation for the density per- 
turbation that can be solved analytically or numerically. 
This yields the transfer function, which determines how 
an initial density distribution evolves as a function of 
time in the presence of small perturbations. 



1. Newtonian theory of small fluctuations 



A further substitution can be made by observing that 
Sp = {dp/dp)6p = cl5p where Cs = (9p/ap)adiabatic is 
the speed of sound. We can also eliminate terms by ob- 
serving that the original (unperturbed) state is spatially 
homogeneous, hence Vp = Vp = 0: 

^ + V • V6p + 6pV • V + pV • (5v = 0, (67a) 
at 



dSv 

— + (v • V)6w + {5w ■ V)v = — '-V5p + (5g, (67b) 
at p 



In order to see how this theory can be developed for 
MOG, we must first review how the density perturbation 
equation is derived in the Newtonian case. Our treatment 
follows closely the approach presented by Weinberg [l^ . 
We begin with three equations: the continuity equation, 
the Euler equation, and the Poisson equation. 



dp 

di 



V • (pv) = 0, 



— + (v-V)v = --Vp + g, 



(64a) 



(64b) 



V • (5g = -AirGdp. (67c) 
Now we note that v = Hx, hence 

V • V = i/V • X 3i?, 

(^v • V)v = {6w ■ V)(i/x) = H{Sv ■ V)x = H6v. 
Therefore, 
dSp 



dt 



V • VSp + 3HSp + p\7 -Sw^O, (68a) 



V • g = -AnGp. 



(64c) 



— + (v • V)(5v + H6v = — '-V6p + 6g, (68b) 
ot p 



First, we perturb p, p, v and g. Spelled out in full, we 



get: 



djp + Sp) 
dt 



+ V ■[{p + Sp){^^ + ^v)]^0, (65a) 
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+ [(v-f(5v) • V](v-H(5v) 



V{p + dp) + g + (5g, 



p + 5p 



(65b) 



V • (5g = -AnGSp. 



(68c) 



The next step is a change of spatial coordinates to co- 
ordinates comoving with the Hubble flow: 



This means 



and 



X = a(t)q. 



V- (g + (5g) = ~4TrG{p + Sp). 



(65c) 



Vq = aV, 
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After this change of coordinates, our system of equations 
becomes 



^ + mSp + -pV • 5v = 0, 

ot a 



— + H6v ^ — '-\/6p + Sg, 
ot ap 



V ■ Sg= -AnaGSp. 



(69a) 

(69b) 
(69c) 



Now is the time to introduce the fractional ampHtude 
S — Sp/p. Dividing (|69ap with p, we get 



S+-S + 3Hd+-W • (5v = 0. 
p a 



(70) 



However, since p = poOg/a'^, and hence p/p = — 3a/a, 
the second and third terms cancel, to give 



a6 = V(5v. 



(71) 



Taking the gradient of (|69bp and using (|69cp to express 
V • Sg, we get 

^(-a6)+H(^aS) = -^V^S-'inGapS. (72) 
ot a 

SpeUing out the derivatives, and dividing both sides with 
a, we obtain 



S + 2H5 - -^V'^S - AirGpS = 0. 



(73) 



For every Fourier mode 6 = Skit)e^^''^ (such that V^(5 
—k^5), this gives 



(5k + 2H6i, 



AttGp Uk = 0. (74) 



The quantity k/a is cahed the co-moving wave number. 

If k is large, solutions to (17^ are dominated by an os- 
cillatory term; for small fc, a growth term predominates. 

A solution to (I74p tells us how a power spectrum 
evolves over time, as a function of the wave number; it 
does not specify the initial power spectrum. For this rea- 
son, solutions to ([71]) are typically written in the form of 
a transfer function 



Tik) 



6k{z ^ 0)^o(z = oo) 
6kiz = oo)(5o(2 = 0) ' 



(75) 



If the initial power spectrum and the transfer function 
are known, the power spectrum at a later time can be 
calculated (without accounting for small effects) as 



P{k)=T\k)P„{k). 



(76) 



P{k) is a dimensioned quantity. It is possible to form 
the dimensionless power spectrum 



A^(fc) = Afc^T^(fc)Po(A:) 



(77) 



where A is a normalization constant determined by ob- 
servation. This form often appears in the literature. In 
the present work, however, we are using P{k) instead of 
A(fc). _ ^ 

The initial power spectrum is believed to be a scale 
invariant power spectrum: 

Po{k) cx fc", (78) 

where n ~ 1. The best current estimate on n is n = 



2. Analytical approximation 

Eq. ([74|) is not difhcult to solve in principle. The solu- 
tion can be written as the sum of oscillatory and growing 
terms. The usual physical interpretation is that when 
pressure is sufficient to counteract gravitational attrac- 
tion, this mechanism prevents the growth of density fluc- 
tuations, and their energy is dissipated instead in the 
form of sound waves. When the pressure is low, however, 
the growth term dominates and fluctuations grow. Put 
into the context of an expanding universe, one can con- 
clude that in the early stages, when the universe was hot 
and dense, the oscillatory term had to dominate. Later, 
the growth term took over, the perturbation spectrum 
"froze" , affected only by uniform growth afterwards. 

In practice, several issues complicate the problem. 
First, the early universe cannot be modeled by matter 
alone; it contained a mix of matter and radiation (and, 
possibly, neutrinos and cold dark matter.) To correctly 
describe this case even using the linear perturbation the- 
ory outlined in the previous sections, one needs to resort 
to a system of coupled differential equations describing 
the different mediums. Second, if the perturbations are 
sufficiently strong, linear theory may no longer be valid. 
Third, other nonlinear effects, including Silk-damping 
[20l |. cannot be excluded as their contribution is signif- 
icant (indeed. Silk damping at higher wave numbers is 
one of the reasons why a baryon-only cosmological model 
based on Einstein's theory of gravity fails to account for 
the matter power spectrum.) 

Eisenstein & Hu [2lj addressed all these issues when 
they developed a semi-analytical solution to the baryon 
transfer function. This solution reportedly yields good 
results in the full range of < fib < 1. Furthermore, un- 
like other approximations and numerical software codes, 
this approach keeps the essential physics transparent, al- 
lowing us to adapt the formulation to the MOG case. 

Eisenstein & Hu [21] write the transfer function as the 
sum of a baryonic term Tf, and a cold dark matter term 

Tel 

Tik)^^nik) + ^TM, (79) 

where fl^ represents the cold dark matter content of the 
universe relative to the critical density. As we are investi- 
gating a cosmology with no cold dark matter, we ignore 
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Tc- The baryonic part of the transfer function departs 
from the cold dark matter case on scales comparable to, 
or smaller than, the sound horizon. Consequently, the 
baryonic transfer function is written as 



Tb{k) 
with 

where 
and 



fo(fc,l,l) a6e-('=/'=siik)^ 



l + (fcs/5.2)2 l + {Pb/ksf 
In (e + 1.8/3,g) 



sin ks 



To{k,ac,Pc) 



In (e + 1.8/3c9) + Cg2 



C 



14.2 



386 



ac 1 + 69.9g 



1.08 ' 



q = kQiT{n^h^y 

The sound horizon is calculated as 



G , \/l + Rd + ^/ Rd + .Rcq 

■ m 



^R. 



(80) 



(81) 



(82) 



(83) 



(84) 



cq 



The scale at the equalization epoch is calculated as 



7.46 X ^0-'''^lr,^h^e2^^ 



(85) 



The transition from a radiation-dominated to a matter- 
dominated era happens at the redshift 



Zeq = 2500017^/1^92^^7, 



(86) 



while the drag era is defined as 



Zd 1291- 



[l + blin„^h^)'^-], (87) 



1 + 0.659(r!™/l2)0.828 I 

where 

bi = 0.313(f7™/i2)~o-4i9[i + 0.607(f]™/i2)"'574]^ (88) 
and 

62 = 0.238(f7„;i2)0.223^ (gg) 

The baryon-to-photon density ratio at a given redshift is 
calculated as 



1000 



(90) 



The Silk damping scale is obtained using 



(91) 

The coefficients in the second term of the baryonic trans- 
fer function are written as 



a, = 2.07fceqs(l + Rdr'^'F ( 

1 + 



(92) 



P{k) 




FIG. 5: The matter power spectrum. Three models are com- 
pared against five data sets (see text): ACDM (dashed blue 
line, ^6 = 0.035, fie = 0.245, Oa = 0.72, H = 71 km/s/Mpc), 
a baryon-only model (dotted green line, Qb = 0.035, H — 
71 km/s/Mpc), and MOG (solid red line, a = 19, p = 
5h Mpc"^ = 0.035, H ^ 71 km/s/Mpc.) Data points 
are colored light blue (SDSS 2006), gold (SDSS 2004), pink 
(2dF), light green (UKST), and dark blue (CfA). 



/3b = 0.5 + + 3 - 2-i^ ) v/(17.2a„/i2)2 + 1, (93) 



where we used the function 
Fiy) = y 



+ (2 + 3y) In 



yr+y+i 



(94) 



A shifting of nodes in the baryonic transfer function is 
accounted for by the quantity 



m = - 



1 + (Aiodc/fcs)" 



1/3- 



(95) 



where 



2\0.435 



(96) 



The symbol 62.7 — T/2.7 is the temperature of the 
CMB relative to 2.7 K, while h = H/{100 km/s/Mpc). 
The wave number k is in units of Mpc~^. 



B. Density fluctuations in Modified Gravity 



We consider the MOG Poisson equation (j48|) . estab- 
lished in section Hi Dl As the initial unperturbed distri- 
bution is assumed to be homogeneous, p is not a function 
of r and can be taken outside the integral sign: 



$y(r) = Gjvap 



1 



-p|r-r'|^3j./^ 



(97) 



10 



Varying p, we get 

V • 5g(r) = -^ttGnSp{y) 

J |r-r'| 
Accordingly, ([75)1 now reads 



2H5 - -4^S/^S - AttGnpS 



(98) 



P^GnchpS 



cPt' = 0. 



(99) 



The integral can be readily calculated. Assuming that 
|r — r'l runs from to the comoving wavelength a/k, we 
get 



7r/2 2tt a/k 



'3„' 



r - r' 



-(fv' 



sin 9 dr d(h dd 





47r [1 - (1 + fia/k)e-t"'^''] 



(100) 



Substituting into ([99)) . we get 



5 + 2i?(5 - ^V^J - AttGnpS 



AttGnch 



or 



1_ fl + ^)e-'^aA- 



(5 + 2HS 



p6 ^ 0, (101) 



-s V72 



V^(5 (102) 



= 0. 



This demonstrates how the effective gravitational con- 
stant 

Geff = G'Ar{l + «[l- (l + ^)e-^'^/'=]} (103) 

depends on the wave number. 

Using Goff, we can express the perturbation equation 

as 



2HSu 



- 47rGeffp Uk = 0. (104) 



As the wave number k appears only in the source term 



AttGcsP , 



it is easy to see that any solution of (fM]) is also a solution 
of ()104|) . provided that k is replaced by k' in accordance 
with the following prescription: 



fc'2 = fc2^4™2/'G«ff-G^ \ ^-2 



Gat 



A7^ (105) 



where Aj = ^/c^/GnP is the Jeans wavelength. 

This shifting of the wave number applies to the growth 
term of the baryonic transfer function ([80]) . However, 
as the sound horizon scale is not affected by changes in 
the effective gravitational constant, terms containing ks 
must remain unchanged. Furthermore, the Silk damping 
scale must also change as a result of changing gravity; as 
demonstrated by Padmanabhan [20], fcsiik c>c G'^/^. Us- 
ing these considerations, we obtain the modified baryonic 
transfer function 



sinfcS j fo(fc',l,l) 
kS ^ j 1 -I- (/cs/5.2)2^ 

g^exp (-[fc/fcsiik]^-VG3/4) 
1 + (/3f,/fcs)3 



(106) 



The effects of these changes can be summed up as fol- 
lows. At low values of k, the transfer function is sup- 
pressed. At high values of fc, where the transfer function 
is usually suppressed by Silk damping, the effect of this 
suppression is reduced. The combined result is that the 
tilt of the transfer function changes, such that its peaks 
are now approximately in agreement with data points, as 
seen in Figure O 

Data points shown in this figure come from several 
sources. First and foremost, the two data releases of the 
Sloan Digital Sky Survey (SDSS) [13, [13] are presented. 
Additionally, data from the 2dF Galaxy Redshift Sur- 
vey (24i], UKST t25!], and CfAl30 |26!] surveys are shown. 
Apart from normalization issues, the data from these sur- 
veys are consistent in the range of 0.01 h Mpc~^ < k < 
0.5 h Mpc~^. Some surveys provide data points outside 
this range, but they are not in agreement with each other. 



C. Discussion 

As a result of the combined effects of dampened struc- 
ture growth at low values of k and reduced Silk damping 
at high values of fc, the slope of the MOG transfer func- 
tion differs significantly from the slope of the baryonic 
transfer function, and matches closely the observed val- 
ues of the matter power spectrum. On the other hand, 
the predictions of MOG and ACDM cosmology differ in 
fundamental ways. 

First, MOG predicts oscillations in the power spec- 
trum, which are not smoothed out by dark matter. How- 
ever, the finite size of samples and the associated window 
functions used to produce power spectra are likely mask- 
ing any such oscillations in presently available survey 
data (see Figure (5]). These oscillations may be detectable 
in future galaxy surveys that utilize a large enough num- 
ber of galaxies, and sufficiently narrow window functions 
in order to be sensitive to such fluctuations. 

Second, MOG predicts a dampened power spectrum 
at both high and low values of k relative to ACDM. Ob- 
servations at sufficiently high values of k may not be 
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1 X Iff*- 




5x 1(P -I • • • 1 • • — 

0.05 0,1 0,15 0,2 0,25 0,3 

k 

FIG. 6: The eflFect of window functions on the power spectrum 
is demonstrated by applying the SDSS luminous red galaxy 
survey window functions to the MOG prediction. Baryonic 
oscillations are greatly dampened in the resulting curve (solid 
red line), yielding excellent agreement with the data after 
normalization. A normalized linear ACDM estimate is also 
shown (thin blue line) for comparison. 

practical, as we are entering sub-galactic length scales. 
Low values of k are a different matter: as accurate three- 
dimensional information becomes available on ever more 
distant galaxies, power spectrum observations are likely 
to be extended in this direction. 

In the present work, we made no attempt to account 
for the possibility of a non-zero neutrino mass, and its 
effects on the power spectrum. Given the uncertainties 
in the semi-analytical approximations that we utilized, 
such an attempt would not have been very fruitful. Fu- 
ture numerical work, however, must take into account the 
possibility of a non-negligible contribution of neutrinos to 
the matter density. 

VI. CONCLUSIONS 

In this paper, we demonstrated how Modified Grav- 
ity can account for key cosmological observations using 
a minimum number of free parameters. Although MOG 
permits the running of its coupling and scaling constants 
with time and space, we made very little use of this fact. 



Throughout these calculations, we used consistently the 
value of a ~ 19 for the MOG coupling constant, con- 
sistent with a flat universe with 0^ ~ 0.05 visible mat- 
ter content, no dark matter, nor Einstein's cosmological 
constant. In nearly all cosmological calculations, we set 
the MOG scaling constant /i to the inverse of the radius 
of the visible universe, which is a natural choice. The 
only exception is the mass power spectrum calculation, 
where the scaling constant enters in conjunction with the 
wave number fc, and describes gravitational interactions 
between nearby concentrations of matter, not on the cos- 
mological scale. 

The theory requires no other parameters to obtain the 
remarkable fits to data that have been demonstrated 
here. 

At all times, lim G — Gn, i.e., the effective gravita- 

r— >0 

tional constant at short distances remains Newton's con- 
stant of gravitation. For this reason, the predictions of 
MOG are not contradicting our knowledge of the pro- 
cesses of the initial nucleosynthesis, taking place at red- 
shifts of z ~ 10^", since the interactions take place over 
distance scales that are much shorter than the horizon 
scale. 

Our calculations relied on analytical approximations. 
This is dictated by necessity, not preference. We rec- 
ognize that numerical methods, including high-accuracy 
solutions of coupled systems of differential equations, as 
in CMBFAST pjj, or TV-body simulations, can provide su- 
perior results, and may indeed help either to confirm or 
to falsify the results presented here. Nevertheless, our 
present work demonstrates that at the very least, MOG 
provides a worthy alternative to ACDM cosmology. 
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